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Abstract 

Kicked double-well system is investigated both analytically and numerically. Phe- 
nomenological formula for ground quasienergy splitting is obtained using resonances over- 
lap criterion in the framework of chaotic instanton approach. Results of numerical calcula- 
tions of quasienergy spectrum are in good agreement with the phenomenological formula. 

Introduction 

Semiclassical properties of systems with mixed classical dynamics is a reach rapidly developing 
field of research. One of interesting results obtained in this direction is a chaos assisted tun- 
neling. It was shown that the structure of the classical phase space of Hamiltonian systems 
can influence such purely quantum processes as the tunneling [HE]. It was demonstrated in 
numerical simulations that existence of chaotic motion region in the classical phase space of the 
system can increase or decrease tunneling rate by several orders of magnitude [T[ |3] . Typically 
one considers tunneling between KAM-tori embedded into the "chaotic sea". The region of 
chaotic motion affects tunneling rate because compared to direct tunneling between tori it is 
easier for the system to penetrate primarily into the chaotic region, to travel then along some 
classically allowed path and to tunnel finally to another KAM-torus [H [5]. 

Chaos assisted tunneling phenomenon as well as the closely related coherent destruction of 
tunneling were experimentally observed in a number of real physical systems. The observation 
of the chaos assisted tunneling between whispering gallery-type modes of microwave cavity 
having the form of the annular billiard was reported in the Ref . [6] . The study of the dynamical 
tunneling in the samples of cold cesium atoms placed in an amplitude-modulated standing wave 
of light provided evidences for chaos-assisted (three-state) tunneling as well [7j . Recently, the 
coherent destruction of tunneling was visualized in the system of two coupled periodically 
curved optical waveguides [8]. 

The most popular methods which are used to investigate the chaos assisted tunneling are 
numerical methods based on Floquet theory [9j. Among other approaches to chaos-assisted 
tunneling we would like to mention the path integral approach for billiard systems ^2] and 
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quantum mechanical amplitudes in complex configuration space |TT]. In this paper we will 
consider the original approach based on instanton technique, which was proposed in [121 [131 E] 
and numerically tested in [TSl [IS] • 

Instanton is a universal term to describe quantum transition between two topologically 
distinct vacuum states of quantum system. In classical theory the system can not penetrate 
potential or dynamical barrier, but in quantum theory such transitions may occur due to 
tunneling effect. It is known that tunneling processes can be described semiclassically using 
path integrals in Euclidean (imaginary) time. In this case instantons are soliton-like solutions 
of the Euclidean equations of motion with a finite action. For example, in Euclidean Yang-Mills 
theory distinct vacuum states are the states of different Chern-Symons classes, and instanton 
solutions emerge due to topologically nontrivial boundary conditions at infinity. In this paper we 
will use a much simpler system to investigate the connection between the tunneling, instantons 
and chaos, namely the kicked system with double well potential. Our attention will be focused 
on the behavior of the quasienergy spectrum when perturbation of the kicked type is added to 
the system. 



1 Instantons in kicked double-well potential 

Hamiltonian of the particle in the double-well potential can be written in the following form: 



Ho = - haox^-aax^ (1) 

2m 



where m - mass of the particle, ao, 02 - parameters of the potential. 
We consider the perturbation 



+00 



Vp,r = ex Sit-nT), (2) 

n=— 00 

where e - value of the perturbation, T - period of the perturbation, t - time. 
Full Hamiltonian of the system is the following: 

H = Ho + Vper- (3) 



Euclidean equations of motion of the particle in the double-well potential have a solution 
- instanton. In phase space of nonperturbed system instanton solution lies on the separatrix. 
Perturbation destroys the separatrix forming stochastic layer. In this layer a number of chaotic 
instantons appears. Chaotic instanton can be written in the following form: 

chaos -^inst ~l~ ^ chaos- 

It is a solution of the Euclidean equations of motion. Here Xchaos and Xinst - chaotic and 
nonperturbed instanton solutions, respectively, Axchaos - stochastic correction. 

It is convenient to work in the action-angle variables. Using standart technique [17] we 
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obtain expressions for this variables in the following form: 



J{E) = ^Cr^ + (Hx) - Jl-^^E K{x) ) , (4) 




2 \ V '^2 



K-\x), (5) 



2ao X 



where J and - action and angle variables, E - energy of the particle, K{x) and L{x) - full 
elliptic integrals of the first and second kinds, respectively, x) - elliptic integral of the first 
kind, where is integral argument, x - modulus which is in the following way expresses through 
energy E 

X 



1 + Jl-AS&E 



We will use in our analytical estimations expression for the action of the chaotic instanton. 
We assume this action is equal to the nonperturbed instanton action corresponding to the 
Euclidean energy which is less than one of the nonperturbed instanton solution. We expand 
expression (jl]) in powers of the (Euclidean) energy difference from the separatrix ^ = Eg^p — E, 
where E - energy, Egep - energy on separatrix, and neglect terms higher than linear one. As a 
result we has the following linear expression 

frri 

Sixchaosir,^)] = nJ{Esep - = S[xinst{r,0)] - a J— ^, (6) 

V ^2 

where 5'[xj„st(T, 0)] = a^"^ / {?> a^) - nonperturbed instanton action, a = (1 + 181n2)/6 - 
numerical coefficient. 

When stochastic layer in Euclidean phase space is formed under the action of the perturba- 
tion a number of chaotic instanton solutions appear, see fig. [U (b). These solutions implement 
the saddle points of Euclidean action functional and thus give significant contribution to the 
tunneling amplitude. In the some sense their role is analogous to the role of the multiinstanton 
configurations in the case of the nonperturbed system. In the framework of the chaotic in- 
stanton approach the width of the stochastic layer determines the contribution of these chaotic 
solutions of Euchdean equations of motion to the tunneling amplitude. The stochastic layer 
width estimated using resonances' overlap criterion. For this purpose we should expand the ex- 
pression for the perturbation ([2D in the series of exponents. At first we evaluate the coordinate 
X of the particle using expression ^ in the following way: 



and rewrite it as a series of sines 



X 



2V2.l^uY, f r_i ^^^[(2ri-lM, (7) 

y QjC\ 1 — Qj 
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Figure 1: (a) Phase space portrait for the kicked double- well potential in the real time, (b) 
Phase space in the imaginary (Euclidean) time with inset where chaotic layer is showed. Chaotic 
instantons are placed in this layer. In both cases closed curves correspond to regular trajectories, 
scattered points to chaotic ones. The model parameters are m = l,ao = 1/128, 02 = 1/4, 
e = 0.01, z/ = 0.5. 

where a = exp[—'7iK{x')/K{x)] - parameter (x' = 1 — u - the frequency of the particle 
oscillations. Leaving only the first term in expression ([7]) and performing simplifications of this 
relation we obtain the following form for the coordinate 



Using this relation and performing the Fourier transformation of the series of delta functions in 
the perturbation one can obtain the following equations of motion in action-angle variables 



where B - angle variable for the particle in the system, = 27r/T is a perturbation frequency. 

To obtain width of stochastic layer we calculated the parameter of the resonances' overlap. 
It is equal to relation of the resonance width in the frequency scale (Acj) and the distance 
between resonances (Su). First one is estimated as a frequency of oscillations of the resonance 
angle variable \E' = — nur. The result is obtained from the equations ([H]) and ([9]) using 
standart technique (see chapter 1 and 5 in [18]). The resonance width is the following 



where AH = Esep — E - distance from the separatrix. 

The distance between resonances is calculating using the expression for the resonance levels 
uJn = n u. Thus one can obtain 





8 



(9) 




(10) 



600 = 00. 



'n-l-l — — 



U. 



(11) 
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Using two last expressions ( ITOl) and ( ITTi) the parameter of the resonances' overlap can be 
written in the following form 



OUJ 



1/2 



> 



e V 
AH 



1/2 



(12) 



Overlap parameter is equal to unity on the boundary of the stochastic layer. Using equation 
(1121) we can write the expression for the width of the stochastic layer in the following way 

AHs = Esep - Ebor ^keu, (13) 

where Ebor is the energy on the border between stochastic and regular regions, k - some nu- 
merical parameter which can not be obtained in the framework of the criterion used. 

The tunneling amplitude for the perturbed system is a sum of the amplitude in the nonper- 
turbed case and the amplitude of tunneling via chaotic instantons. The later can be evaluated 
by integration over action of the tunneling amplitude in nonperturbed system. Using expression 
([6]) this integral can be transformed to the integral over the energy difference from zero up to 
the width of the stochastic layer (fT3l) : 




Achaos = 0!\—N / dco\/ S[xchaos{r, 0] ea;]? {-S[xchaos{r, C}]), 



where is a normalize factor. To calculate contribution of chaotic instantons we use approxi- 
mate expression for the chaotic instanton action ^ and listed above assumptions. Integration 
over Co gives the contribution of zero modes |19J. As the result we get the following expression 
for the amplitude: 



A = Ainst + A^haos ~ N e"^'"'* TexplaJ— AH^ ) , (14) 




where Amst is tunneling amplitude in the nonperturbed system, F - a time of the tunneling which 
is put to infinity at the end. The last exponential factor in the expression (H^ is responsible for 
the tunneling enhancement in the perturbed system. In the nonperturbed case the width of the 
stochastic layer is equal to zero and the expression (fT^ coincides with the known expression 
describing the ordinary tunneling. 

Substituting AHs in formula ( JT^ by the expression (fT3l) and expanding it in powers of the 
perturbation strength we can write phenomenological formula for the quasienergy splitting 

Ar7(e,i/) = 2y?v/5^e-^'""(l + A;ei/), (15) 

where 

k = Oa — k. 
V (^2 

We fix this phenomenological parameter value using the results of numerical simulations. For 
this purpose we perform the linear fitting of the numerical data for the dependencies on the 
perturbation strength and take average value of the parameter over these dependencies. As the 
result we have the single numerical parameter k for our phenomenological formula explaining 
all these dependencies. 
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2 Numerical calculations 



For the computational purposes it is convenient to choose as basis vectors the eigenvectors 
of harmonic oscillator. In this representation matrix elements of Hamiltonian ([T]) and the 
perturbation ([2]) are real and symmetric. They have the following forms (n > m): 



m+4 n 



aog 



A/(m + l)(m + 2)(m + 3)(m + 4) 



+5 

m n 



+ 5m+2 n^ig ao(2m + 3) - 4) A/(m+ l)(m + 2) 
1 /3 

huj{n + -) + I ( - ^ ao (2m2 + 2m + 1) - a'2(2m + 1) 



e 5. 



m+l n 



vm+T, 



where (7 = h/muj and a'g = 02 + mco'^/2, /i - Planck constant which we put equal to 1, a; - 
frequency of harmonic oscillator which is arbitrary, and so may be adjusted to optimize the 
computation. We used the value u = 0.2 with parameters m = 1, oq = 1/128, 02 = 1/4. In 
numerical simulations size of matrices was chosen to be equal to 200 x 200. Simulations with 
larger matrices give the same results. 
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Figure 2: Quasienergy spectrum for the ten levels with the lowest average energy. The model 
parameters are m = 1, ao = 1/128, 02 = 1/4 and u = 0.5. Thick lines - doublet with the 
minimal energy. 

We calculate eigenvalues of the one period evolution operator e~*^"^e~*^ and obtain quasienergy 
levels which are related with the evolution operator eigenvalues through the expression 77^ = 
i InAfc/T. Then we get the ten levels with the lowest one period average energy which is calcu- 
lated using the formula {vi \ Hq + V/T \vi) {\vi) are the eigenvectors of the one period evolution 
operator). The dependence of quasienergies of this ten levels on the strength of the perturbation 
is shown in the figure [21 Quasienergies of two levels with the minimal average energy (thick lines 
in the figure [2]) has a linear dependence on the strength of the perturbation in the considered 
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Figure 3: Quasienergy splitting as a function of the strength (a) and frequency (6) of the per- 
turbation. Lines - phenomenological formula, points - numerical results. The model parameters 
are m = 1, = 1/128, a2 = 1/4. 

region. They are strongly influenced by the perturbation while some of the quasienergy states 
are not. 

Performed numerical calculations give the dependence of the quasienergy splitting on the 
strength (figl^a)) and the frequency (figj3]^b)) of the perturbation. These dependencies are 
linear as it predicted by the obtained analytical formula (fT^ . Using least square technique we 
calculate the single numerical parameter k which describes all numerical dependencies demon- 
strated in the figures [3]( a) and[3](b). Relative error in determining of the parameter k from 
numerical results is less than 0.1%. Analytical results are plotted in the figures [3] (a) and (b) 
by lines. Numerical points lie close to this lines. The agreement between numerical simulations 
and analytical expression is good in the parametric region considered. 



Conclusions 

Double-well system is investigated in presence of external kick perturbation. Analytic chaotic 
instanton approach is applied for this system in order to obtain the phenomenological formula 
for the ground quasienergy splitting. The formula has the single numerical parameter which is 
determined from the numerical results. This formula describes ground quasienergy splitting as a 
function of strength and frequency of perturbation. It predicts linear dependence of the ground 
quasienergy splitting on these parameters. Numerical results for the quasienergy splitting as 
a function of the perturbation frequency and strength demonstrate linear dependence as well. 
They are in a good agreement with the analytical formula. Analytical determination of the 
parameter k at the moment seems to be impossible due to estimative character of the analytical 
method used. However, it is not a large problem since this parameter can be evaluated with the 
high accuracy using the single dependence of the ground quasienergy splitting for any particular 
system parameters' values. Parameter k determined in such a way allows to describe correctly 
all other dependencies of the ground quasienergy splitting on the perturbation strength and 
frequency for any other values of the system parameters. The only restriction is that the 
perturbation strength has to be sufficiently small (in the case considered e < 0.1) in order the 
chaotic instanton approach used to be valid. 
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